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A complementary code pair has the property that its autocorrelation 
function contains no range sidelobes. This property makes complementary 
codes attractive for use in signal waveform design for detection and 
tracking systems. An extensive list of the properties of these codes 
has been formulated by M J. Golay and additions were made by S. Jauregui. 
However, the ambiguity functions of complementary codes have not been 
previously formulated • 

This thesis considers the ambiguity functions of complemientary 
codes for two cases: 

1. Narrow-band analysis 

2, Wideband analysis. 

In the narrow-band analysis it is assumed that the doppler effect 
does not significantly alter the envelope of the transmitted signal. 
Siebert*s definition of the ambiguity function is utilized in the formu- 
lation of the narrow-band complementary code ambiguity function. 

In the wideband analysis it is shown that the doppler effect sig- 
nificantly alters the envelope of the transmitted signal. This alteration 
of the envelope of the transmitted signal becomes significant when the 
product of the two factors, the relative velocity between the transmitter 
and the target, and the time-bandwidth product, approaches one-half the 
velocity of propagation in the medium, Siebert's ambiguity function is 
extended to include the wideband case. 
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CHAPTER I 



INTRODUCTION 

The detection problem in modern search radar systems is generally 
concerned with the optimization of four major factors. These fouf 
factors ares 

1) Maximization of the effective range of the system, 

2) Minimization of the number of false alarms due to 
impulse noise o 

3) Maximization of range and velocity resolution, 

A) Accurate measurement of the target’s position and velocity 
without ambiguity. 

The solution of this optimization proglem is predicated primarily 
on a suitable choice of a transmitted signal waveform. To emphasize 
this point let us analyze the four factors that are involved in the 
optimization problem. 

The first factor can be analyzed from the standpoint of the radar 
range equation. The prominent result of this analysis is that the 
system designer may increase the range of the system by increasing the 
duration of the transmitted pulse. 

The second factor can be disposed of by setting a high enough 
threshold such that very few noise impulses exceed the threshold, but 
not so high as to degrade the probability of detection below a tolerable 
level. 

The last two factors depend almost entirely on a judicious choice 
of a signal waveform. Hence it is topical to consider a measure or cost 
function by which various waveforms may be compared, Siebert’s^ modifi- 
cation of Woodward's ambiguity function provides such a comparison. 
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It is for this reason that the ambiguity function forras the basis of 
the modern radar detection philosophy* 

This thesis is primarily concerned v;ith a class of codes, called 
complementary codes, x^jhich minimizes range ambiguities* !Iost codes, 
which have been described and analyzed in the literature on radar wave- 
form design, have the drawback that range and doppler ambiguities exist. 
Examples of these codes are: 

1) Linear RI (chirp) 

2) Truncated pseudo-noise codes 

3) Barker codes. 

Of the codes mentioned above, the Barker codes* autocorrelation functions 
exhibit the smallest range sidelobes, and have been named "perfect*' words. 

Complementary codes have the feature that their autocorrelation 
functions have no range sidelobes* Although much work has been done in 
the area of obtaining properties of these codes, the analysis of the 
ambiguity functions for these codes has not been done. It is the pur- 
pose of this thesis to analyze complementary codes, in both range and 
doppler, for tx>7o cases: 

1) Narrow-band analysis 

2) Wideband analysis* 

By narrow-band we mean that the doppler effect does not significantly 
alter the envelope of the transmitted signal. The narrow-band analysis 
will be formulated by utilizing Siebert*s ambiguity function.* lioxcever, 

*Siebert*s ambiguity function is a modification of Woodx<^ard*s 
ambiguity function* The essential difference is that Siebert's ambiguity 
function is real valued, whereas Woodward’s ambiguity function is complex 
valued . 
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one of the postulates in the formulation of Siebert’s ambiguity function 
is that doppler does not significantly affect the envelope of the trans- 
mitted signal. Hence Siebert*s ambiguity function cannot be used in the 
x^ideband analysis. We xv^ill extend Siebert*s ambiguity function to include 
the wideband case® In Appendix B we x^7ill define the region of validity 
of the narrox7-band analysis. The essence of the derivation of the region 
of validity of the narrow-band analysis is that the timie-bandwidth-velocity 
product must be mucli less than the velocity of propagation divided by 
tv7o, in order that we may assume that the envelope of the transmitted 
signal is not altered by the doppler effect. Then by x^7ideband, x^7e mean 
that doppler has a significant effect on the envelope of the transmitted 
signal. Throughout this thesis, \^e will assume that target acceleration 
is negligible over the interval of measurement. 

Let us now consider complementary codes. Complementary codes or 
series were originally conceived by M. J. Golay,^>^ in his x^ork on 
infrared multi-slit spectrometry. An extensive list of properties for 
complementary codes was formulated by H. J. Golay and additions, based 
upon applications of group and vector theory, x^ere formulated by 

c 

S. Jauregui in his doctoral dissertation.*^ 

A set of complementary series may be defined as *'a pair of equally 
long, finite sequences of two kinds of elements which has the property 
that the number of pairs of like elements x-/ith any one given separation 
in one series is equal to the number of pairs of unlike elements x^7ith 
the same given separation in the other series.**^ These sequences are 
binary and are aperiodic. 

It has been shown that the autocorrelation function of these codes 
gives rise to a central peak. The height of the peak is dependent on 
the length of the code, for zero delay, and decreases linearly to zero 
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for delay equal to the duration of one code syrbol, then remains at zero* 
It is this property of the conplementary codes that makes this coding 
method attractive for use in the detection problem* Because the defi- 
nition of the complementary codes appears rather formidable, v;e will 
utilize the property given above whenever we refer to complementary 
codes* Appendix F is a partial list of the properties of complementary 
codes. This list is provided so that complementary codes of almost any 
length may be generated from basic length codes called kernels, A kernel 
is defined as a basic length code which cannot be decomposed into shorter 
length codes by an inversion of the standard generating methods. To the 
author’s best knowledge there are only four kernels that are known to 
exist 0^ These kernels are listed in Appendix F. 

At this point we might indicate what the ideal ambiguity function 
would be. Since we would like to resolve two targets in range and doppler 
to an infinite degree, the ideal ambiguity function would be a delta 
function at the origin of the range, doppler axis. Thusly, the system 
would have infinite resolution. This can never be achieved; what would 
be required is a code of infinite length with the duration of each code 
symbol infinitesimally sm.all. However, we should not discard the idea 
of an ideal ambiguity function, since it gives us a reference by which 
to compare other ambiguity functions. 

A brief synopsis of each chapter will be now given. 

In Chapter II a system model for single codes is postulated and a 
correlation between a physical waveform in the system and the ambiguity 
function is given. By single codes we are referring to binary codes 
such as Barker codes, pseudo-noise codes, etc., i.e., codes which con- 
sist of a single sequence; as opposed to complementary codes which 
consist of a pair of equally long binary sequences. In this chapter 
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the analysis is entirely narrow-band. The single code aniLipuity function is 
derived to ease computational matters in the derivation of the ambiguity func- 
tion for complementary codes. 

Chapter III contains the derivation and the justification of the form of 
the narrow-band complementary code ambiguity function. 

In Chapter IV a system model is postulated for the complementary codes 
and the derivation of the complementary code ambiguity function is accomplished 
Computer programs were run for the two kernels of length 10 and the kernel of 
length 26, and ambiguity diagrams were constructed. One property of the entire 
class of binary codes in which phase-reversal modulation is utilized, is un- 
covered in this chapter. Tt is proved that in the doppler domain, the zero 
crossings of the ambiguity function move closer to the origin of the range- 
doppler axis as the length of the code is increased; assuming that the duration 
T, of a code symbol remains constant. The implication of this phenomenon is 
that better doppler discrimination is possible for longer code lengths, all 
other things being equal. 

Chapter V contains the investigation of the doppler effect on the trans- 
mitted signal and the proof that the envelope of the transmitted signal is 
altered due to doppler. This chapter in conjunction with Appendices A and D 
shows that the doppler effect is threefold: 

1) The carrier frequency is shifted by the doppler frequency. 

2) The envelope of the transmitted signal is compressed or expanded 
according to whether the target is closing or opening with respect 
to the receiver. 

3) The received and transmitted energies are different; in fact, the 
received and transmitted energies differ by the same factor as the 
received and transmitted carrier frequencies. The above statement 
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assumes a lossless medium and ignores spherical spreading of 
the transmitted and received waves. 

Chapter VI contains the derivation of the wide-band ambiguity 
function for single codes as well as complementary codes, valid for 
positive differential doppler and positive range deviation. Unfortunately 
due to time limitations, computer programs to generate numerical values 
for the wideband complementary code ambiguity function were not accom- 
plished. 

To clarify the discussion of the model of the system, which is 
postulated in Chapter II, a brief formulation of the autocorrelation 
function of complementary codes will be given. As mentioned previously, 
complementary codes consist of two binary seqtiences, which will be 
referred to as Code A and Code B or simply A and B. To obtain the auto- 
correlation function^ for the complementary code, we compute 

the autocorrelation function for code A, and for code B, Rg(r). 

Then R^(T ) is obtained by summing ) and i.e., 

The property of interest in the autocorrelation function of the comple- 
mentary codes is that no range sidelobes exist. All binary codes have 
autocorrelation functions with range sidelobes. However, has 

range sidelobes which are the negative of the sidelobes of R^(T). 

Hence when R^(T^) ^^d are summed, the range sidelobes cancel. Fig- 

ure 1-1 shows this unique property of the complementary codes. 
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Fig. 'i-ia. - Autocorrelation function for code A, length 10 
A = set ^-1 ,^1 ,+1 ,-1 ,+1 ,-1 ,+1 ,+1 ,+1 ,-1} 

Fig. 1-1b - Autocorrelation function for code B, length 10 
E - £.t i --' .+1 , + 1 , M ,+1 ,+1 ,-M ,-1 ,-1 ,+ 1 } 

Fig. i-'c - Aui ocori-eiati- fanction for the complementary 
r*^ae, l‘‘rgth 10. 
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CHAPTER II 



MATCHED FILTER APPROACH FOR THE 
SINGLE CODE SYSTE-l (NARROW- BAND ANALYSIS) 



It might be well to look at the single code system before attempting 
to analyze the complementary coded system. In the matched filter 
approach the underlying theme is to correlate the ambiguity function 
with a physical waveform in the system. The ambiguity function is an 
excellent measure of range and velocity resolvability, i,e., how closely 
one can distinguish two targets in range and velocity. A comon mis- 
conception is to state that the output of the matched filter is identical 
to the ambiguity function for all cases. This statement holds some 
truth; however, it is not genera[|.ly correct in carrier m.odulated systems. 
Let us depart from this discussion to model our system before attempting 

to clarify the above statement. The model that is chosen is a slight 

8 

modification of the model used by Davenport. 



Transmitter 






J Transmitted 








Medium 



s * (t-t d) -f n(t) Receiver 



'4 

information 
source 
(target) 



I disturbance I 
(additive) 



System model 
Figure 2-1 

The difference between Davenport’s model and the model that is chosen 
is that in the radar, as well as the sonar case, the information source 
interacts in the medium, whereas in Davenport’s system model, the infor- 
mation source interacts in the transmitter. 
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The transmitted signal is chosen as a phase-reversal modulated 



signal, where 



where 



s(t) = \cos(u"Qt - n^" ) 2-1 

^i ” J ^ depending on the value of the ith code 
' 1 

symbol 

\ = normalization constant . 



For convenience the signal energy is normalized to unity by specifying 
that 

U^(t) dt = 1 

It is a well known fact that a phase-reversal modulated signal is 
identical to a DSB/SC signal; hence the transmitted signal can be written 
as 

s(t) = I(t) cos ^u^t, 2-2 

r ^ H 

where I(t) = V (t-iT) - _it - (i + l)Tlj^ 

and = - 1, depending on the value of the ith code 

symbol* 



The normalization constant can be easily evaluated by squaring equation 
2-2 and performing the integration* Thus 




(t) dt = 




NT 

\ 1 + cos2 Lv' ^t ] dt 

J o ^ o 



-2 9 

■X NT rj (^sin2 iOqNT 

2 4 'Co 



2-3a 



2-3b 



The second term in equation 2-3b is very small compared to the first 
term, for large* This term can be interpreted as the difference 
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between the initial and final value of the carrier phase and will be 



small. Hence the normalization constant is approximately given by 



0 ( = 






2 

NT 4 



2-3c 



The received signal contains a noise term, n(t), which will be 

n 

ignored since the receiver will utilize a matched filter to maximize 
the signal-to-noise ratio. Thus the received signal, s’(t-t^), is the 
transmitted signal modified by the information source, and contains 
range and velocity information. Let us assume that the target speeds 
are sufficiently less than the speed of propagation of the wave in the 
medium (say of the order of 10~^ x c, v/here c is the velocity of propa- 
gation,) The principal doppler effect, in this case, is to shift the 
carrier frequency by the amount f^, where f^ is the doppler frequency; 
doppler will have little effect on the duration of a code symbol, T. 

(See Chapter IV for the discussion of the doppler effect on the duration 
of a code symbol.) With the above conditions postulated, the received 
signal may be given by 

cos[Ci4,+ 



j'j 



where the convention is chosen so that the doppler frequencv is positive 
for closing targets, and negative for opening targets. 

At this point the question might be raised as to v;hat are the 
various v/aveforms in the receiver. As mentioned earlier in this ch«apter, 
a matched filter will be utilized in the system m.odel. To obtain the 
signals to correlate with the received signal, v;e translate the trans- 
mitted signal by various expected doppler frequencies as indicated in 
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Figure 2-2a 




XCr.uwj 

4 

tc Comparator*^ 



'^Envelope 

detector’’ 



Figure 2-2b 

Figure 2-2a Receiver block diagram 

2-2h ith doppjer filter of the filter bark 
in Figure 2-2a 
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Figure 2-2a. The doppler shifted signal is tlien delayed. The delay is 
required to account for the propagation time to and from the target. To 
instrument this system in real time^ the delay mechanism might be of the 
form of a tapped delay line. We shall denote the delayed, doppler shifted 
signal as the replica signal. The replica signal is given by 

Si(t - = I(t - ric)cos [(Oio + - Zyj] 2-5 

Let us now analyze the various signals in the receiver. To simplify 
the analysis, the matched filter block diagram is duplicated from 
Figure 2“2b, 

Tlie output of the matched filter is 

-oo 






Then let 



= t - tj 



>^01 " *^0 + ^di 



Equation 2-6b becomes 



2-7a 

2-7b 






^ ±(gi(i^4 Cos(u}e,-vO^. 

“•OO 



Cos[ Woj.(-t^4-ta-'rk^l , 2-8 
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Then let 




2.9a 



(Jld “ h)(ji 



2.9b 



If the above changes in variables are effected, then 



OO 




2-10 



The variables T and deserve some comment. From equation 2-9a, 
we see that if is zero, then the propagation delay to the target 

Tk 

has been matched in range bin Tt . Hence — is the range to the 

^ 2 

target. From equation 2-9b we see that if (jC^d zero, then Q)di is 
the doppler frequency. Thus both Z snd represent the departure 
of the target’s range and velocity from the estimated range and velocity. 
The variable t^ is just the variable of integration, and since it is 
inconsequential as to what nomenclature is ascribed to it, we shall 
denote by t. cOq! represents the estimated carrier frequency of the 
replica signal. Hence the output of the matched filter in the system 
is given by 




Kt'iKt- r) cos u)oj.(t-r)dt ^ 



2-lla 



-00 



60 




00 




2-llb 
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The second term in equation 2-llb has frequencies centered about twice 
the carrier frequency plus the doppler frequency, and will be rejected. 
Note that the output of each matched filter is an ordinate in the ^ 
domain. 

Next we require a device to collect the output of each of the 
matched filters in an orderly fashion. Hence the use of the term, 
collator, to describe that device is appropriate. We shall define the 
collator in Figure 2-2b as that device which collects the various out- 
puts of the matched filter sequentially as a function of range deviation, 
. To clarify the above definition, we might think of the collator as 
a graph plotter, where ordinates (outputs of the matched filters) are 
plotted as a function of range deviation, T • The “envelope detector'* 
in Figure 2-2b is not an envelope detector in the normal sense, and 
requires some explanation. As mentioned above, the collator collects 
the outputs of each matched filter and produces a continuous curve in 
the X, domain. This curve oscillates at in the X domain, but only 

its envelope has information content. The “envelope detector" refers 
to that device that takes the envelope of the output of the collator, 
but in the range deviation or T domain. In the usual sense of the 
term, envelope detection is accomplished in the time domain and hence 
a distinction must be made in our usage of the term envelope detector. 

Now let us look at the output of the collator; 



If the Euler identity for the cosine function is utilized, the output 
from the collator is 




2-12 
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1 

= IT® 









4 i ^ 2U)ia-r)e^‘^‘^^dt 



2-13 



If we note that the second tern of equation 2-13 is the conplex conjug:ate 
of the first term and that the sum of a quantity, say z, and its complex 
conjugate is 

Z + = 2 Ile(Z) ^ 

then 



= Z j/e {i 



2-14 



The terms I(t) and I(t - T ) are the modulation functions and involve 
the limits of integration and the values of the codes. The modulation 
functions, according to equation 2-2 are 



M 

= 0( >1 X; I iCt-iT)- ii (■ 2-15a 

0(2 Xj^ilVZ-lT)- 2-15b 



The unit step functions in equations 2-15a and 2-15b A^7ill deterir.ine the 
limits of integration and the X^'s are ll, depending on the value of the 
codes. Figure 2-4 is the plot of the time waveforms of I(t) and I(t - r ). 
To obtain the limits of integration, merely observe the zero crossings 
of the functions. Hence the output of the collator is given by 
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Figure 2-4 
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The above equation can be generalized to a code of any length I’, and 
will be written as the sum of two series. Therefore the output of the 
collator is 
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where N 
k 
T 



length of the code 
smallest integer 2 1 t/np \ 

the duration of a code symbol = clock period 



The first integral in equation 2-16 can be evaluated as 
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2-17b 



” ai^lK-vilT + j siu uoiCk+OT - Gosaj^Cr-i-Ci-DT] 

-j SiD. (i-i)Tl "*1^ , 

If the following trigonometric substitutions are utilized, 
cos X- co£ y = - 2. sia sin 
siax ± sii\Y = '2 ^ 



equation 2-17b becomes 
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The second integral in equation 2-16 can be evaluated as 
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By similar arguments used to derive equation 2-18d, becomes 
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We then substitute equations 2-18b and 2-19b into equation 2-16, 
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If the exponential terms are combined, and the real part of equation 
2-20 taken, then 
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In order to plot equation 2-21, it is necessary to obtain the equation 
in the form of an envelope term multiplied by an RF terra. Let 
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We now utilized the trigonometric identity 



cos (a + b) = cos a cos b - sin a sin b , 



Equation 2-23a becomes 
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Thus the output of the collator becomes 
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where 0 is the phase due to the amplitude of the sine and cosine terms. 



If equations 2-22 a, b, c, d, e, are substituted into equation 
2-23c, and if the sine terms are placed in the familiar form^ 
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then 
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Equation 2-24 is plotted in Figure 2-5; note the RF fine structure 
under the envelope. 

Now let us analyze the ambiguity function. Essentially, definitions 
and results will be stated from Siebert’s paper^ with the exception of 
a few terms which will be modified to conform to the nomenclature of 
this thesis. The ambiguity function as defined by Siebert is 
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where S^(t) is the complex envelope of the transmitted signal. 
In terms of this thesis 

S(t) = I(t). 



I(t) is real. Therefore 
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Equation 2-26 departs from equation 2-14 in that the carrier phase term 
is absent and the magnitude of the integral is taken. Utilizing the 
same analysis used in obtaining the limits of integration fpr the output 
of the collator, we see that equation 2-26 can be written as 
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figure 2-5 

Waveform of the output of the collator for the 
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If the results of equations 2-17 through 2-20 are utilized, equation 
2-27 becomes 
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The term e 






can be factored from equation 2-28. Note that 



e ^ 



= i 



If the magnitude of equation 2-28 is taken, then the equation in its 
final form is 
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2-29 
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In the earlier portions of this chapter | the statement was made 
that the output of the matched filter is identical to the ambiguity 
function. We are now prepared to discuss this particular point. In 
comparing equation 2-24 (the output of the collator) with equation 2-29 
(the ambiguity function), we note that the ambiguity function is identical 
to the envelope of the output of the collator. Hence the difference 
between the ambiguity function and the output of the matched filter is 
that there exists an RF fine structure of frequency f i iri the output 
of the matched filter whicli is absent in the ambiguity function. Note 
that if the phase of the RF fine structure is some odd multiple of 
then the output of the matched filter will be zero. Thus, we can say 
qualitatively that the ambiguity function measures the positive raxinuns 
of the output of the matched filter per cycle of the RF fine structure. 

To circumvent the problem of outputs of the matched filter being zero 
near the central peak, we sum and take tlie envelope of the output of 
the matched filters. However, care must be taken in choosing the differ- 
ential range deviations. Suppose that the differential range deviation, 
/iT f is some multiple of tl\e inverse of twice the frequency of the RF 
fine structure, i.e., 
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Furthermore, suppose that in any matched filter, the phase difference 

between the transmitted si[,nal and the received signal is an odd multiple 

of JL- radians 0 Then the output from all the matched filters x^/ill be 
2 

zero. This point is easy to see because the system will be tracking the 
zero crossings of the RF fine structure. 

The RF fine structure is also troublesome in the area of multi- 
plexing. This point will be discussed in Chapter IV, when we will 
formulate ideas on the type of multiplexing to be used in the complemen- 
tary coded system. 

The output of the "envelope detector" is given by the envelope of 
equation 2-24. However, as mentioned earlier, the output of the "enve- 
lope detector" is identical to the ambiguity function. 
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CHAPTER III 

DEFINITION OF THE AMBIGUITY FUNCTION 
FOR COMl^LEMENTARY CODES 

In this short chapter, we shall extend Siebert*s Arbiguity function 
to the ambiguity function for compleirentary codes. Let us analyze the 
signal processing techniques required to produce the desired results in 
the complementary coded system. In Chapter II, it was shown that in 
the single code system, the output of the ^'envelope detector” in Figure 
2-2b was identical to the ambiguity function. Uliat then, should be the 
configuration of the complementary coded system? Should the complementary 
coded system be, in essence, two single coded systems such that the 
outputs of the two "envelope detectors” be summed together? Or should 
the two codes be summed elsewhere? 

Let us postulate the complementary coded system to be com.prised of 
two identical single code systems with the outputs of the "envelope 
dectors” of code A and B summed together. In this case the ambiguity 
function of the complementary coded system will be the sun of the ambi- 
guity functions of code A and B, i.e., 

One of the properties of the ambiguity function is that 






= 0 ( 
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Hence equation 3-1 for UD = 0 becomes 



1^=0 



3-2 



As indicated in the introduction of this thesis, the desired property 
of the complementary codes was that off the central peak, the autocorrelation 
function of code A, was the negative of the autocorrelation function 

of code B, so that when the autocorrelation functions were sumjned, 

the range sidelobes canceled « Since the proposed definition (equation 
3-1) of the ambiguity function contains two positive quantities • 

and definition is inadequate. 

As a second approach to the problem, let us define the ambiguity 
function of the complementary coded system as 





3-3 



where = envelope of code A 

IgCt) = envelope of code B 



For zero doppler frequency, equation 3-3 becomes 




3-4 
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It is easy to see that the range sidelobes of code A V7ill cancel the 
range sidelobes of code B so that the range sidelobes for the comple- 
mentary coded ambiguity function will be nonexistent for zero doppler 
frequency. Equation 3-3 will be chosen as the definition of the ambig- 
uity function for the coirpleirentary codes. 

Equation 3-3 gives us the clue as to the configuration of the 
complementary coded system. We can infer that the two codes should be 
summed at the output of the matched filters, then envelope detected. 
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CHAPTER IV 



MATCHED FILTER APPROACH FOR THE COMPLEMENTARy 
CODED SYSTI?! (NARROW- BAND ANALYSIS) 



A . BACKGROUND 

In the complementary coded system, some form of multiplexing is 
required since two codes will be utilized to obtain the desired cancel- 
lation of the range sidelobes® Of the three obvious types of multi- 
plexing (time division, frequency division, and quaternary coding) only 
time division, TDM, and frequency division multiplexing, FDM, will be 
discussed. 

B. PROPOSED FDM SYSTEM 

Assume for the moment that the FDM system is modelled as below. 
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Complementary coded FDM 
transmitter block diagram 

Figure 4-1 

It might be well to observe the transmitted spectrum to aid in the 

understanding of the signal processing schemes that v/ill be required in the 
receiver. The frequencies, and f 2 » are at RF, The oiatputs of each 
of the phase-reversal modulators can be easily obtained from equation 2-2 
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if appropriate subscripts are used. Hence 



tJX) = 0(^ 2 XA-, \ila- iT) - i[l - (U 1)T] CCS , 

1=1 

iT 

0(,2X6,(ik.t-irt-lIU-(i^l)T] , 

® i*l ^ 



It is easy to see that the spectrum of the spectrum of the transmitted 

sin X 



signal consists of two 



functions centered at and f 2 > 




One-sided spectrum of the transmitted signal 
Figure 4-2 



Thus to detect each code, the system requires bandpass filters centered 
at the carrier frequency plus the doppler frequency. Hence the system 
block diagram is given by Figures 4-3a and 4-3b. Since the output of 
each range bank can be treated the same as in the single code system, 
equation 2-15 applies with the substitution of appropriate subscripts. 
Let be the envelope quantity in equation 2-15; hence the output from 
any range bank in Figure 4-3b for code A is given by 
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FDM receiver block diagram 
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4-3a 



The output from the corresponding range bank for code B is 




4-3b 



The output from the summer in Figure 4-3a is the sum of the outputs 
from the range banks of code A and B, Therefore 



The effect of the different frequencies of the RF fine structure for 
code A and B is to lower the output from the matched filter at the 
central peak and to increase the range sidelpbes. It is true that, if 
the carrier frequencies of both codes are nearly equal, the range side- 
lobes will be very small; and in fact, might be quite tolerable* However, 
hear in mind that the characteristic of interest in a complementary 
coded system is for range sidelobes to be nonexistent for zero doppler* 

In this system the factor which is troublesome is the RF fine 
structure t^hich is absent in the ambiguity function. Hence the analysis 
provides a further contradiction to the statement that the output of the 
matched filter is identical to the ambiguity function. Indeed, it is 
that factor which is different between the output of the matched filter 
and the ambiguity function which defeats this method of multiplexing. 
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C. TDM SYSTEM 

Another multiplexing method is TDM. The suceeding portions of 
this thesis will dwell entirely on phase-reversnl TDM systems for comple- 
mentary codes. The results from Chapter II concerning single code 
systems can be easily extended to the complementary coded systems. In 
the transmitter the essential difference is the time multiplexing circuitry; 
hence the transmitter can be modelled as below. 




cos (^^t 



TDM transmitter block diagram 
Figure A-4 

As in Chapter II the transmitted energy will be normalized to unitv 
by specifying that 
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To obtain the normalization constant, we note that 
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s^(t) and Sj^(t) are given by 



- W = XA;{iU-iT) - l[l- (Ui)tlcos a)A , 

Uo ^ ^ " 
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where XA^ = Zl^ depending on the value of the ith code symbol of A 
XB^ = tl, depending on the value of the ith code symbol of E* 
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Therefore the normalization constant, O^c* approximately 



0^. 
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In the transmitter, code B is delayed and summed with code A, The 
receiver is modelled in Figure A-5, To detect the signal the gating 
circuit in the receiver separates the two codes* The output from each 
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TDiV; receiver block diagram 
Figure 4-5 
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correlator in Figure 4-5 is identical to the output from correlator 
in the single code system. Hence the output from the ^envelope defector** 



is given by 



= Envelope of {y^(T^ -v y^Cr, ^ 
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Upon examination of equation 2-15, we find that the only terms which 
depend upon the value of the codes arc the The rest of the terms 

are common to both y^ ( T , Vb ( T , ^d^* if the signed 

products in equation 2-15 are replaced by signed products for code A and 

B, i.e., 



and if the envelope of equation 2-15 is taken> then the output of the 
**envelope detector'* will be given by 
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As stated in Chapter II the output from the "envelope detector" is iden- 
tical to the ambiguity function. Therefore the ambiguity function for 
the complementary codes will be given by equation 4-9. 

Let us now analyze the narrow-band ambiguity function for comple- 
mentary codes. For zero doppler equation 4-9 reduces to 



Thus we see that from the analysis done in the introduction, the ambiguity 
function for zero doppler contains no range sidelobes. We note that 
this property is a direct result of the RF phase synchronization assumed 
in equations 4-6a, b. 
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For the case of zero range delay, the ambiguity function reduces to 
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Note that 



X.X. = 1 . 
1 1 



Then if the sine terms are factored, equation 4-11 reduces to 
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If the Euler identity is utilized 
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and noting that 
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then, we see that 
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Thus equation 4-12 can be written as 
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If the following trigonometric forraulas^^ are utilized 
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thus equation 4-13 becomes 



a. 






z 



KuiAT QOS t;jg 0 

a z 



,CK+i^(AlT 

s -m — 



Then noting that the absolute value of a product is the product of 
absolute values , and 
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Finally we see that equation 4-14 becomes 
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If equation 4-15 is plotted, one interesting result is observed 




Ambiguity function for complementary codes^ 
in the doppler domain. 

Figure 4-6 

The first zero crossing of the ambiguity function occurs at succeeding 

zero crossings occur at multiples of Hence the coded pulse sequence 

N1 

is acting like a single pulse of length NT in the doppler domain. For 
a better insight into this phenomena, let us review the origin of the 
sin term. A glance at equation 4-12a shows us that the term origi- 

nated as a sum of exponentials, i.e.. 
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The criterion for the ambiguity function to be zero in the doppler domain 
is 



_ Q 
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Or alternatively 



ja = A, 2^. ^ 4-17b 

Hence we see from equation 4-16 that when the total phase shift is a 
multiple of 2T ^ the ambiguity function for zero range delay will be 
zero. This effect can be viewed as the addition of phasors such that 
when the function is zero, the phasors have added to form a closed loop. 
This, of course, occurs when the total phase shift is some multiple of 
2tt , 

One added comment should be made concerning the zero crossings of 
the ambiguity function in the doppler domain. Since the zero crossings 
occur at then, as the length of the code increases, the ambiguity 

function will be compressed closer to the origin; hence approaching the 
ideal. However, the length of the code cannot be increased indefinitely 
without complications# It will be shown in Chapter VI that longer codes 
will be more sensitive to changes caused by the doppler effect, and 
unfortunately the length of the code cannot be increased indefinitely 
without reducing the output from the matched filter. The above compression 
effect on the ambiguity function is not a property possessed solely by 
the complementary codes. For example, single phase-reversal codes have 
ambiguity functions which experience the compression effect for long 
codes. It is the length of these codes, rather than their internal 
structure, which produces compression in the doppler domain; assuming 
that the duration of a code symbol, T, is held constant. 
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D. SUMMARY OF RESULTS 



Equation 4-9 was coded into Fortran and programs were run for the 
two kernels* of length 10 and the kernel of length 26. Plots of the 
ambiguity function were made for constant *Moppler phase^'* ^ where 

this is defined as the recurrent argument in our equations for the ambi- 
guity function, i.e.. 



A 

where T = duration of a code symbol. 

To convert doppler phase into radial velocity, we recall that for 
the narrow-band case 





Therefore^ 



A A 

2* io 



C 

^ ioTrr , 
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Table 1 gives the conversion from doppler phase to radial velocity for 
various values of carrier frequency and duration of a code symbol. 
Doppler phase, , was defined in equation 4-18^ so that the ambiguity 

function need not be recomputed for different values of f^ and T. 



*As stated in the introduction, a kernel is a basic length code 
which cannot be decomposed into shorter length codes by an inversion of 
the standard generating methods. Conversely, complementary codes of 
longer lengths can be generated by appropriately combining the kernels. 
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Figure 4-7 is the ambiguity diagram for the kernel, and B 2 , 

code length 10 ^ (The code values for the kernels are listed at the end 
of Appendix F.) The ambiguity diagram is the term used to describe the 
three dimensional figure generated from the ambiguity function. The 
bandwidth* of the ambiguity function in Figure 4-7 is approximately 
O.lTT . For the system postulated in Figure 2-2, this bandwidth corresponds 
to a velocity discrimination of 3000 knots, for T = 1 ju^sec, f^ = 10 GHz. 

For f^ = 1 GHz and T = O.l^sec., velocity discrimination is roughly 
300,000 knots. Hence complementary codes of this length cannot be used 
for situations requiring good velocity discrimination. 

Figure 4-8 is the ambiguity diagram for the kernel, A 2 and B 2 , 
code length 10, and with in increments of 0.125'U' . 

Figure 4-9 is the ambiguity diagram for the kernel, and B^^, 

code length 10. Note that the range sidelobes extend out to seven units, 
whereas the range sidelobes in Figure 4-8 extend to 10 units. However, 
the amplitudes of the sidelobes for the kernel, Aj^ and Bq^ are slightly 
larger than for sequences A 2 and ^ 2 ^ 

Figure 4-10 is the ambiguity diagram for the kernel, A and B, length 
26. The range sidelobes extend to 14 units. The maximum amplitude of 
the range sidelobes is 0.32. The ambiguity diagrams in Figure 4-9 and 
4-10 are in the same increments of . Note that the bandwidth of 

the ambiguity function for the kernel, A and B, length 26, is less than 
the bandwidth of the ambiguity function for the kernel, Aj^ and length 
10. This effect is the compression phenomena for longer codes which is 
discussed in Chapter IV. 



*The bandwidth of the ambiguity function is here defined as doppler 
phase (converted to frequency) at the first zero crossing of the ambiguity 
function for = 0. 
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TABLE 1 



CONVERSION FROM DOPPLER PHASE, » TO 

TARGET’S RADIAL COMPONENT OF VELOCITY, 

Note: c = velocity of propagation of the wave in the medium 

(5*91 X 10° knots) 

f^ = carrier frequency 

T = duration of a code symbol 

f^ = doppler frequency 

= TTfjT = 2vrfoTTr/C 



v^ (knots) 


fg (OIIZ) 


T ( sec , ) 


^ j (tt rad) 
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TABLE 1 
(cont 'd) 



(knots) 


fg (GHz) j 


T ( jusec.) 


(TTrad) 


60 


10 


1 
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CHAPTER V 



INVESTIGATION OF THE DOPPLER EFFECT 
ON THE TRANSMITTED SIGNAL (WIDEBAND ANALYSIS) 



A. BACKGROUND 

In Chapters II, III, and IV the analysis was done on pulse-coded 
systems for the narrow-band case* l\Jhat is meant by narrow-band is that 
doppler has negligible effect on the envelope of the transmitted signal, 
but shifts the carrier frequency by the doppler frequency. It is the 
purpose of this chapter to discuss the doppler effect on the envelope 
of the waveform and to derive the expression for the time waveform of 
the received signal (wideband case). As stated in the introduction, 
target acceleration is assumed to be negligible over the interval of 
measurement* 

B. EFFECT OF DOPPLER ON THE TRANSMITTED SIGNAL'S FREQUENCY COMPONENTS 

Assume that s(t) is any arbitrary transmitted signal. Then the 
Fourier transform (if it exists) is a measure of the frequency content 
of the signal. It is assumed that s(t) is sufficiently well behaved 
to h'^vo a Fourier transform. If the shorthand notation is utilized, the 
Fourier transform of s(t) is 



It is well known that the doppler effect is to shift the frequency 
components of the signal such that the new frequency f' is given by 



3" Lscnl = . 




5-1 



*This equation is a statement of the two-way doppler effect for 
radar and sonar. For the radar case Torres^l has indicated that if the 



where = radial component of velocity between the target and the 
transmitter 

c = speed of propagation of the \<rave in the medium 
Note: The convention is chosen so that v^ is positive for closing 

targets and negative for opening targets. 

Define a quantity^ the doppler factor^ such that 

f* = Df ^ 5-2a 



Alternatively, 



0 5-2b 

A-Vf/c • 

The doppler factor may be given in terms of the doppler frequency f^. 
Note that the doppler frequency is defined as the change in the carrier 
frequency due to the doppler effect. Hence 

fd = fo' - ^o . 5-3b 

Then 

^d = - fo = (D - l)f„ . 5-3c 



velocity vectors of the transmitter and the target are noncollinear , the 
angle at which the energy arrives at the receiver (assuming that the same 
antenna is used for transmitting and receiving) differs from the angle 
at which the energy was transmitted due to the angular aberration phe- 
nomenon.^^ 

In the sonar case equation 5-1 applies if the medium is at rest 
with respect to the transmitter. 
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The above equation nay be rearranged so that 



•D = 1 + 



5-4a 



If the ratio — « 1, 



tlien the doppler factor may be approximated by 



D = ± + 

c 



5-4b 



Hence the doppler frequency is approximately 



f 



d 



2vr 

c 



5-5 



Next let us examine the received signal’s Fourier transform, S’(f). 

(The prime notation is used to indicate that the received signal differs 
from the transmitted signal®) Since the effect of doppler is to shift 
all the frequency components of the original signal by a multiplicative 
factor, D, the received signal’s transform is given by 

S'(f) = S(f/D) 5-6 

C. TIME WAVEFORM OF THE DOPPLER SHIFTED SIGNAL 

It was established in the previous section that the received signal’s 
Fourier transform was related to the transmitted signal’s transform by 
a change in argument. Hence the inverse transform of the received sig- 
nal’s Fourier transform will yield the time waveform of the received signal. 



*We arrived at the amplitude factor being unity from considerations 
of the special theory of relativity. In Appendix A we proved that the 
transmitted and received energy is dependent on the relative velocity 
between the transmitter and the target. This result was also obtained 
by A. W. Rihaczek. To obtain the energy of the received signal, we 

integrate the squared magnitude of equation 5-6 over all frequencies. 
Thus the energy amplitude factor is introduced via the change in fre- 
quency scale. 
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Since the time waveform and its Fourier transform form a one to one 
correspondence, the time waveform obtained by the inverse transform of 
the received signal's Fourier transform v;ill be unique. If the inverse 
Fourier transform is taken, then 



-oO 



By a change in variable 




then equation 5-7 becomes 



oo 



m = o\ 



5-8 



Equation 5-8 has the same form as the transmitted signal with the excep- 
tion of the amplitude factor D and the change in the time variable to 
Dt. Therefore 



Ds(Dt)- 



Equation 5-9 is a well known transform pair and is referred to by 
13 

R. M. Bracewell as the "Similarity Theorem." Since the received sig- 
nal is of the same functional form as the transmitted signal except for 
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a change in argument and in amplitude, s'(t) may be written as 

S'(l^ = D J(04») COS 5-lOa 

= T)0<2^Xi{A(D^-'li-il'LDi-(i+l)T]}cos 5-lOb 

Note that the unit step function jL(Dt-lT) may be written as and 



Define 






'o= '/d. 



5-11 



Then 






5-12 



Hence equation 5-lOb becomes 






s'(-t) = Oo< { la-lTo') . 

i = i 



5-13 



Equation 5-13 is the time waveform for the wideband doppler shifted 
signal. Note that the wideband signal reduces to the narrow-band signal 
if the doppler factor is nearly unity. 
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If equation 5-13 is examined carefullyi we see that doppler produces 
three effects, 

1) The carrier frequency is shifted by the doppler frequency, 

2) The duration of a code symbol is altered by the doppler factor, 

3) The received energy differs from the transmitted energy by the 
doppler factor. 

The reason for the carrier frequency being shifted by doppler is obvious, 
and Appendix A shows the relationship between the received and transmitted 
energy. However, why is the duration of a code symbol altered by the 
doppler effect? 

Let us determine the physical significance of the alteration of the 
duration of a code symbol. 

Let 

^ = number of cycles of RF in the duration of one code symbol 
of the transmitted signal 

= number of cycles of RF in the duration of one code symbol 
of the received signal. 



Then 



fi = duration of one code symbol for the transmitted signal 
^ duration of one cycle of the transmitted signal 



= ^ = f T 

T7fT 



5-l4a 

5-l4b 



_ duration of one code symbol for the received signal 5-15a 
duration of one cycle of the received signal 



l/fo’ 



= fo’To . 



5-15b 



From equation 5-11 



T 

T = — 

D , 
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Time waveforms of the transmitted and received 
signals for the duration of one code symbol 





69 



and from equation 5-3b 



f • = Df 



o • 



o 



Therefore 



= (Dfo)(n) = foT 



5-16 



Hence 



'I - ^ 



5-17 



Thus we see that the number of cycles of RF in the transmitted and 
received signals remains constant* Hence the implication of the doppler 
effect is that when the transmitted signal's carrier undergoes a shift 
due to doppler, the duration of a code symbol changes so that the total 
number of cycles of RF for the transmitted and received signals remains 
constant* This is compatible with the special theory of relativity as 
shown in Appendix D. 
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CHAPTER VI 



WIDEBAND AMBIGUITY FUNCTION FOR 
COMPLEMENTARY CODES 



A. BACKGROUND 

It might be well to indicate those situations in which the doppler 
effect on the envelope of the transmitted signal cannot be ignored. 
Certainly for relative velocities, between the target and the transmitter^ 
near the velocity of propagation of the wave in the medium, this effect 
cannot be ignored. Another situation involves the duration of the pulse 
train. As indicated in Chapter V, the doppler shifted waveform undergoes 
a compression or expansion; hence if the total duration of the shifted 
waveform is compressed or expanded, approaching the duration of one code 
symbol, energy will be lost in the correlation process. The criterion 
for the use of the narrow-band analysis is given in Appendix B and is 
duplicated here. Hence if 

« ^ 6-1 




Or alternatively, 



NTWlv^l I 



(See equation B-14) 



6-2 



Then the narrow-band analysis applies. 



*This result is indicated by A. W. Rihaczek^^, and E. J. Kelly and 
R. P. Wishner,14 
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Equation 6-2 is a statement of the velocity-time-bandvidth-product 
uncertainty relationship, 

B. SYSTEM BLOCK DIAGRAM 

The system block diagram for the wideband case is essentially the 
same as the narrow-band case with the exception of two critical areas, 

A guard interval is necessary between code A and code B to allow for 
time expansion or compression of the two codes due to the doppler effect 
on the envelope of the received signal. Also the delayed replica of the 
transmitted signal must be clocked at the rate corresponding to the 
expected clock rate of the received signal and the carrier frequency must 
be shifted by the doppler frequency. Suppose that the replica signal 
were shifted by the doppler frequency, but not clocked at the expected 
rate; then degradation of the output of the matched filter would be 
expected. In fact, the above statement implies that we are requiring a 
narrow-band system to detect a wideband signal. From the point of view 
of equipment complexity, it is more economical just to shift the carrier 
frequency of the transmitted signal to obtain the replica signal. Appendix 
C contains the analysis of the replica signal just being shifted in 
carrier frequency for the wideband case. It is expected that no serious 
degradation of the output of the matched filter, and hence the “envelope 
detector,** will be evidenced, so long as the time-bandwidth-product-velocity 
uncertainty relationship is not violated. (Equation 6-2) 

C. WIDEBAND AMBIGUITY FUNCTION FOR SINGLE CODES 

In the wideband case the replica signal is obtained by clocking 
the transmitted signal at the expected rate and translating this signal 
by the expected doppler frequency. For convenience the wideband matched 
filter block diagram is given below. 
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Wideband matched filter 
Figure 6-1 



The replica signal is given by 



Si(Dil - - KDjt-rjj') cos , 



6-3a 



where = expected doppler factor 




fo 

^di = expected doppler radian frequency^ 

The received signal is given by 

= D' CjOS [(u^o+ (Wj'l > 

where D' = target's doppler factor 




fo 



= target's doppler radian frequency. 



6-3b 



6-4a 



6-Ab 
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Then the output from the matched filter is given by 



CO 

‘^di •> ‘Ai') \ d-l 

— OO 

= D'ri[D'(4--ta)lcos[(05,+ uid)(Ud)]- 

-~ao 

%'il Co4K->^diHl- ^)]at _ 6-5a 

Now let us shift the time reference by letting 

ti = t - t^ . 

Thus equation 6-5a becomes 



4 = o'iKo'OIlDil^i+ld- cos(u),4u:)^')-l^ 0 

-OO Vi 

cos[(03,+u3^.^C-t^4l^- ^)] ^ 6-5b 



Then define 






6-6 



Since is the variable of integration, it may be replaced by t with 
no loss in generality o Hence the output of the matched filter may be 
written as 



OO 

vj. = 0'\ i(D'-t)xLOiU-rj^] cos(u}„+uOaH 






6-7a 



74 





The second term in equation 6^7b will be rejected, provided that 

• As in the narrow-band case, the output 
of the matched filter is an ordinate in the 'C domain. Then, the collator 
collects these ordinates in the HT domain. The output of the collator 
is given by 



If the Euler identity for the cosine function is used, equation 6-8 becomes 




6-8 



Let 



Woi = 1^0+ ^c\t 



6-9a 



% , 



6-9b 




6-10 



The output of the "envelope detector" will be given by 
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To the extent that we can ascribe an envelope to the output of the 



summer, then 




6-11 



In Chapter II, it was proved that for the narrow-band system the output 
of the '’envelope detector" is identical to the ambiguity function. For 
the wideband case the ambiguity function will be defined by 



It is only to the extent that the approximation made in equation 6-11 
holds, that it is possible to give a physical interpretation to the 
wideband ambiguity function* 

Let us now reference the received signal to the replica signal by 
letting 



*If the integral in equation 6-11 is a "relatively" slowly varying 
function of , as compared with » then this approximation 

is valid. 






6-12 



= Dit. 



Thus equation 6-12 becomes 




6-13 



76 



Then define 



_ o' _ 




6-lA 


Oi 




= - 




6-15 


- _ 
Di 




6-16 


= Co^- 




6-16b 



Therefore the wideband ambiguity function becomes 







,60 

-QO • 
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Note that either equation 6-17 or equation 6-13 may be chosen as the 
definition of the wideband ambiguity function. 

One of the properties of the narrow-band ambiguity function is that 
it is symmetrical in range deviation, T , and doppler, However, 

the wideband single code ambiguity function seems to exhibit no simple 
symmetry property. Thus, we will obtain equations for the wideband 
single code ambiguity function, explicit for a limited case^ to indicate 
the method that is involved. 

The wideband single code ambiguity function v/ill be calculated for 
T>0. These equations will be limited to the case of the 
received signal’s envelope being compressed or expanded by not more than 
the duration of one code symbol of the replica signal. Appendix E 
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contains the analysis of the limits of integration and the code values 
in the interval of integration for the wideband ^ single code amoiguity 
function® Ke have divided the ambiguity function into two intervals of 
range deviation. 

1. For range deviations, ^ , in the interval 



AT 



AT at 






we have 



A 



i-i. 



^ ei' 









6-18 



where N = length of the code 



= il depending on the ith value of the code symbol 
Dj = replica signal’s doppler factor = 1 + 

lo 
fa ’ 

D’ = received signal’s doppler factor = 1 + — ^ 



D = £1 . f p +■ 



D,- + ^di 



V 



6-19 



^A = (D. - 1) Ulo 



-'A 

T = transmitted signal's duration of a code symbol 



"Di 



Td' 



T' = 



replica signal’s duration of a code symbol = — 

Di 

received signal’s duration of a code symbol = ^ 

Ti/ 



78 



AT = mcignitude of the difference in the duration of a code symbol 



= T . - T* 

^Di ^ 



p = smallest integer > 
j = ^-1 (denotes a complex quantity) 
0l = normalization constant ^ 



6-19 
(cont . ) 



Then if the integration is performed, equation 6-18 becomes 












a 



a ^ 



ui 



' A -T3 



a 



C.r^lA^b+lTo.-(i4P^V3 
lr 4 irp.-ius/^T ] ^ — 

i=i ^ o 



a 



‘^A(.rfkTp;-(.i^p>T3 

a 
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2. For range deviations, , in the interval 



£sir 



AT AT AT ’ 



m = 1,2,**,N I 



we have 



^'(r,uia = 



D^0(^ 

a 



M-k LHk)T . 

ZXiXi,K\ 



T+a-DTja 






^r+iTo; 






i^t 
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+ £> dt 

ftW/ 



gfi-Ox ' 



K-2. 

+ 2 

!=«-!! 



T4G-k>T5; 



e> -dt 






(aHK' 






U4Z^t' 

IrO-i j 

rv(i-k^TQ. 



where the definitions in equation 6-19 apply and 

T' 



n = 



AT 



"1 ■ Tt 



Tpi 

A' 



q = smallest integer > 
1 = q modulo (n^^) 
k = smallest integer > 



Vat 






If the integrals in equation 6-21 are evaluated, then the ambiguity 
function becomes 



'I'CTi (AJ = 0^(X^ \ 2 XiXi.i, \1U^)T- (M^Tp;-rl 






^^UUk^T-a-lVToi-Tl 

2x 



i “ — r 



2 . 



l=i 



-21 
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sia 






A 



a 



^ -r'sin'^' 

^ <^aT" 






+ Z Xi.kXi^^ [r4g-kFpi-U4i)T^J 

1-n-C 2^ 



uiA^^+(i-liVoi-.(■a+l^T'Jl . uk[U-kF«4(i4A')T3 

Si JTi *** ^ 

^ACT-t-U-k^~S; - (i-ViKj 



+ X Xi« K»»^)T-(i4JtFi>i-'r] o 

1 = n-1 ^ 



. . (jO^ r(UZ>T - (J4j<>‘^i - t3 

SiJfl 5 



^At(uz)T - (Mk)v.-r] 



. u3AKU2)T4-(i-k^TK3 

i 
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D. WIDEBAND AI4BICUITY FUNCTION FOR COMPLEMENTARY CODES 

As stated in Chapter IV-C, to obtain the ambiguity function for a 
complementary code, given the ambiguity function for single codes, one 
merely replaces the signed product for the single codes by the signed 
products for the complementary codes# 

1. Consider range deviation, T , in the interval 



AT AT At at 



m = 1,2,‘*‘,N. 
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Then if the signed products for the single codes in equation 6-20 are 
replaced by the signed products for the complementary codes , the ambi- 
guity function for the complementary codes will be given by 





o 



z 






O 



e) ^ 



M-Cf+i') 

2 feA;XA; 






1+fA.i 



It4 1 >q. ^ 

£ 



^ 

Z 






> 
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where n/ = 
c 



1 __ 

NT 



2* Now consider range deviation, 'T , in the interval 






Then if the signed products in equation 6-22 are replaced by the signed 
products for the complementary codes, the ambiguity function for the 
complementary codes will be given by 
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= D^CX, 






2 IxAiXAuk 

X^X 



-. ..,„uiAUuk)T-a-lVDl-'r] 
CU ^ toT - (i-DToi-^l z 



L ^^ c ( i + k ') T '- a - i')%-'d 



; » ^ACCUk>T'-+ ( i -J 0 ~^ 0 i 1 

1=1 



-i- XB;X6i,t,.i 1 + iToi - (i-^k'>T'3 « 

Z 



sifi *^ACt^ iTt)L-a^k>T'3 

^Atr 4 lTo|-(l 4 k')T '3 

2 



. L 0 i^ti''cii + (i+W^ i 3 
gj — z 






u^aV ; uWt^Zrv- 2Sl+l)"r- ‘C’3 
^ 

LOaT 

2 

4:2 

”^2a LXA-5_^XAi<.2, X-2>i.i^X6)i+2i"i ° 

i=r,-ii 
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2, 



t ^AQr-vCi-k^'^^-dA 1 )t'J 



. tcu [H-k>Tbi + ( 1 + DT'J xi-2. 

2 



+ 2 






2 



uiAEU2)ri(i4k'r^/-'r] . uD^.LCi-vz'^Tl 
sin-^ 2 — 

^ADi+2>T- Ci-^kFoi-'Tcl 
2 



6-24 
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CHAPTER VII 



CONCLUSIONS 



A • SUM>^RY 

The research presented in this thesis had two primary objectives. 

The first objective was to formulate the ambiguity function of the comple- 
mentary codes for the narrow-band case. The second objective was to 
formulate the ambiguity function of the complementary codes for the 
wideband case. We shall summarize the results of these investigations, 
but not necessarily in the order presented in the thesis. 

1. NARROW- BAND ANALYSIS 

In the narrow-band case, Siebert*s ambiguity function was extended 
to include complementary codes. The ambiguity function of complementary 
codes was defined so that its autocorrelation function contained no range 
sidelobes. From the definition we determined that the two codes should 
be combined at the outputs of the matched filters to effect the necessary 
cancellation of range sidelobes. During the course of research into the 
ambiguity function of complementary codes, we discovered that, once the 
ambiguity function for a single codes has been derived, then to derive 
the ambiguity function for a complementary code is a simple task — replace 
the signed products of the single code x^ith the sum of the signed products 
of sequences A and B, i.e.. 



replace (X^Xj) by C^A^XAj + XB^XBj ) . 

In the derivation of the ambiguity function for single codes, it 
was proved that the ambiguity function was identical to the output of 
the "envelope detector." (The "envelope detector" refers to tlie device 
in the system that takes the envelope of the signal, but in the range 
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or T domain. This use of the term differs from conventional terminology 
in that, usually ^ envelope detection is done in the time domaint) Although 
this result is not unique, it helps to clarify some of the misinterpre- 
tations of the ambiguity function. It is a common nisconcept icr to state 
that the output of the matched filter is identical to the ambiguity 
function. 

One result obtained in the analysis of the ambiguity function of 
complementary codes is that in the doppler domain, the zero crossings 
of the ambiguity function move closer to the origin as the length of the 
code increases, if the duration of a code S3onbol is held constant. The 
implication of this result is that longer codes permit better velocity 
d^‘«^rimination. However, this process of increasing the code length to 
obtain better velocity discrimination cannot be continued indefinitely 
without incurring a penalty. The penalty is that if the velocity-tine- 
bandwidth-product uncertainty relationship is violated, then distortions 
of the envelope of the received signal must be taken into account. This 
amounts to having to clock the replica signal at the expected rate and 
hence increasing the complexity of the system. 

2. WIDEBAND ANALYSIS 

In the wideband analysis, it can be shown that the doppler effect 
is three-fold. 

a. The carrier frequency is shifted by the doppler frequency. 

b. The envelope of the transmitted signal is compressed or 
expanded according to whether the target is closing or 
opening with respect to the receiver. 

c. The received energy differs from the transmitted energy 

by the same factor as the received and transmitted carrier 
frequencies . 
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Although all three effects can be proved from the special theory of 
relativity, only the last two were proved in this manner, The distortion 
of the envelope of the received signal was also proved as a logical con- 
sequence of the application of Fourier analysis, and the doppler effect 
on the carrier frequency. A. W. Rihaczek^^ proved, by considering a 
photon, that the received and transmitted energies differ. Our approach 
to the problem was to make a Lorentz transformation of the field tensor. 

It was proved that for binary phase-reversal modulated signals, if 

NTWlvrl « I , 

then the narrow-band analysis applies. 

The \>7ideband ambiguity function was derived here; hox>7ever this 
result had been previously accomplished by A. W. Rihazek^^^ Kelly 
and Wishner^^. We merely clarified certain aspects of the analysis. 

Then the equations for the single code and complementary code wideband 
ambiguity functions were derived. These equations are limited to the 
case of signal's envelope being altered by not more than the duration 
of one code symbol. They are, indeed, formidable, and a computer is 
required to obtain the ambiguity diagrams. Due to time limitations, the 
wideband ambiguity diagrams were not completed. 

B. AREAS OF POTENTIAL STUDY 

Many areas of future research follow directly from this thesis. 

A logical extension of this thesis would be to construct the wideband 
ambiguity diagrams for complementary codes and single codes (Barker codes, 
truncated pseudo-noise codes, etc.) Also the region of validity for the 
narrow-band analysis should be more precisely defined. We found that, 
in the preliminary analysis of the computer runs for the wideband ambi- 
guity function of complementary codes, the wideband ambiguity function 
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for complementarv codes began to deviate substantially from the narrov- 
band ambiguity function at approximately 

NTwlvrl= 0.1 c/2 
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APPENDIX A 



ENERGY OF THE RECEIVED SIGNAL 
DERIVED FROM CONS IDEPvAT IONS 
OF THE SPECIAL THEORY OF RELATIVITY 



Consider the reference frame, S, with the transmitter located at 
the origin of the coordinate system. Assume that the target is a perfect 
reflector, which is moving at a velocity, v^., relative to the transmitter. 
The target is situated at the origin of the reference system, S*, A 
suitable set of axes is chosen so that the x direction of both of the 
systems, S and S’, is collinear with the relative velocity vector, v^. 
Figure A-1 reflects a suitable choice of axes. Note that in Figure A-1, 
the convention for positive velocity is away from the transmitter and is 
a standard convention for most physics textbooks. Although this con- 
vention is contrary to the one chosen in the main text of this thesis, 
no confusion should arise here. 

The transmitter emits a plane wave of finite duration, N'f, Let field 
quantities as seen in the transmitter reference frame be Ey and H^ , 

Next, we require the field quantities as seen by the target in 
reference system. S’, These quantities can be obtained by application 
of the Lorentz transformations of the field-strength tensor^^. They are 



given by 



= 0 

Ky' = - e>Hz) 

E^' = = 0 



»x' = = 0 

Hy* = 'g'(Hy + ^E^) = 0 A-1 

llj,' - "bCHz - ^Ey) , 



where 



^ = Vr/c 
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s 



S’ 




Reference frames for the target and the transmitter 

Figure A-1 






The plane wave after reflection 
Figure A-2 






The plane wave after reflection 
as viewed by the transmitter 

Figure A-3 
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Then we want the reflected wave’s field components as viewed by 
the target. These field components can be easily obtained by applying 
the boundary conditions for a perfect reflector. Therefore, the field 



quantities as seen by the target are given by 




Hx' ' = 0 




Ky- = 0 



A-4 



E " = 0 



H '• = +11 ' . 
z z 



z 



Finally, we require the field quantities as viewed by the transmitter. 
The reflected v/ave’s components can be easily obtained by applying the 
transformation rules in equation 1, and noting that the velocity component 
is now in the negative x direction; hence we must replace v^ by -v^.. 

Let us denote the reflected field quantities by a triple prime notation. 



Thus 



L ” ’ = 0 

X 




Ly’" = X(Ey" + 
= 0 



1 ! " ' = 0 
y 



A-5 







§ 



If equation A-4 is substituted into equation 





^(-Ey’ + ^ll^') 

^(ll^' - P>Ey')^ 



A- 6 



Then substituting equation A-1 into equation A-6 



Ey’" = '^{[-^(Ey - - f.Iy)} 



= - + 2 

- §Ey)- ^^(Ey - ^U^)] 



A- 7 




A-« 
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The solutions for Maxvv^ell’s equations are such that 



since 




E =112. 



= H, 



= 1 , if the field quantities are in Gaussian units. 



Therefore, equation A-7 becomes 



Ey’" = I- ^^(1 + + 2 

= A-9 

and 

_ 2 ^ 2 ^ 

= ^$^(1 - ^ )2h^ . A-10 

The power density can be obtained by Poynting’s Theorem, and is 
given by 



P=i<5^e(i E.H-i , 

Thus the power density of the reflected wave as seen by the transmitter is 









" L = IS > 



A-12 



v/here is the amplitude of E^. 
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The power density of the transmitted wave is 



Ip I - A- I p 1^ A-13 

The energy density of the wave can be determined by integrating 
the power density over the duration of the signal. Note that the power 
density is uniform over the planar surface. Thus the energy density of 
the transmitted wave is 

^ A-l« 

-oo 



And the energy density of the reflected wave is given by 




A-15 



v;here N ~ length of the code 

T’*’ = duration of a code symbol for the received signal 

T = duration of a code symbol for the transmitted signal* 

During the analysis presented in Chapter V and Appendix D, it is sliown 
that the received signal’s duration of a code symbol is altered by the 
doppler effect so that 

T 1 1 1 2L 

D ' 

where 

D = 
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Thus equation A-15 becomes 



& 



m 








A-16 



Then the ratio of the received signal's energy density to the transmitted 
signal's energy density is given by 



= D . A-17 
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APPENDIX B 



CRITERION FOR THE USE OF THE 
NARROW-BAND APPROXI^^ATION 

(VELOCITY-TIME-BAIT)WIDTH-PRODUCT UNCERTAINTY RELATIONSHIP) 

In this appendix we shall derive, from rather elementary consider- 
ations, the region of validity for the narrow-band approximation* This 
approximation is made in Chapters II and III and assumes that doppler 
has negligible effect on the envelope of the signal. Chapter V contains 
the derivation of the time waveform of the doppler shifted signal for 
the x^ideband case. (For purposes of clarification, the terms, doppler 
shifted signal and the received signal, are used synonymously.) We con- 
cluded tliat the effect of doppler was threefold: 

1. The duration, T^^, of a code symbol for the doppler shifted 
\^7aveform vjas altered so that 

Tj3 = T/D. B-1 

2. The carrier frequency x^7as shifted by the doppler frequency. 

3. The received signal’s energy density differed from the trans- 
mitted signal’s energy density by the doppler factor. 

If the doppler shifted signal’s duration of a code symbol is altered by 
doppler so that the total duration of the signal is expanded or compressed 
by the duration of a code symbol of the transmitted signal, x^e expect 
degradation of the output of the matched filter. In particular, if the 
code is completely random, then for zero range deviation the output of 
the matched filter x^7ill be zero. Hence \^e shall use as our criterion of 
the narrow-band approximation, that the total duration of the received 
signal must be compressed or expanded by much less than the duration of 
a code symbol of the transmitted signal, T. Let us examine two cases. 
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1. Case I: The total duration of the received is compressed by 

much less than the duration of a code symbol of the transmitted signal, 
(This statement implies that positive doppler is present,) 

Figure B-1 indicates the envelopes of the transmitted signal, I(t), 
and the received signal, I(Dt), for some arbitrary code, Pence for the 
condition that the total duration of the received signal is compressed 
by much less than the duration of a code symbol of the transmitted signal, 
the following inequality applies : 

NTj3 » (N - 1)T , E-2 

If the above inequality is rearranged, 

T/Tj3 (N/(N - 1) . E-3 

From equation B-1 

T/Ti3 = D « K/(N - 1) , 

Therefore 

1 + fd/fo - 1) • 

And 

f^/f^ « 1/(K - 1) . E-4 

Alternatively, 



1 + V /c 

1 — « K/(N - 1) . 

1 - Vj./c 



If K » 1, then 



Vj-Zc « 1/2N . 



E-5 
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Time waveforms of the received and transmitted 
signals for the case of the received signal's 
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Note that the tiine-bandx>rid th product is equal to Ih* 

Hence 

KTW 4^ c/2 . B-6 

2, Case II: The total duration of the received signal is expanded 

by much less than the duration of a code symbol of the transmitted signal. 
(This statement implies that negative doppler is present.) 

Figure B-2 indicates the envelopes of the transmitted signal and 
the received signal for some arbitrary code. Hence for the condition 
that the total duration of the received signal is expanded by much less 
than the duration of a code symbol of the transmitted signal, the following 
inequality applies: 

NT » (K - 1)Td . B-7 

The above inequality may be rearranged so that 

T/Tj) » (N - 1) /N . B-8 

Note that the doppler frequency is negative in this case. Hence the 
doppler factor may be written as 




*Ke have defined bandwidth as the video bandwidth to the first null 
of the rectangular pulse spectrum, i.e., W = 1/T. This is approximately 
equal to the ‘‘effective” RF bandwidth. Thus NTW = N. The time-band- 
width product, as used by Torres, was defined in such a fashion (using 
standard deviation in both time and frequency domains) as to be a complete 
and precise description for a Gaussian pulse. If rectangular or other 
shaped pulses are used, a numerical factor of the order of magnitude of 
unity may be introduced into his precise inequality. Hence there exists 
an apparently discrepancy of a factor of 2 between Torres* inequality 
and our inequality B-6. 
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If equation B-9 is substituted into the inequality E-8> then 





f 



1 



B-10 



Alternatively j the doppler factor may be given by 



\vr\ 



D 



c 



1 



B-11 



c 



Hence, substituting equation B-11 into the inequality B-8, then 




E-12 



If the length of the code is much greater than unity 



Then, we note that the time-bandwidth product is equal to the length 
of the code, N# Hence 



Thus if the time-bandwidth product multiplied by the radial component 
of relative velocity is much less than one-half the velocity of the wave 
in the medium, then the narrow-band assumption applies. Inequality B-14 
is a statement of the velocity— time— bandwidth— product uncertainty relation- 
ship. 



NTUT \v^| 4< I . 



B-14 
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APPENDIX C 



EFFECT OF DOPPLEPv MISMATCH, DUE TO THE PvEPLICA 
SIGNAL ONLY BEING TRANSLATED IN FREQUENCY, ON 
THE OUTPUT OF THE ’’ENVELOPE DETECTOR” 



1 . BACKGROUND 

For the wideband case, the system block diagram contained range 
banks, such that in one range bank, the received signal is correlated 
with an exact replica of itself • Note that the replica signal has the 
same carrier frequency and clock period as the received signals Hence, 
in generating this signal from the transmitted signal, the system must 
contain a clocking device, such as a shift register, as well as a fre- 
quency translator • Naturally tlie clocking device adds to the complexity 
and hence the cost of the system* Suppose that in trying to reduce the 
complexity of this system, we merely translate the frequency of the 
transmitted signal by the doppler frequency* How is the output of the 
’’envelope detector” affected? If the replica signal is obtained by a 
mere frequency translation, we are, in effect, requiring a narrow-band 
system to detect a wideband signal* Thus we expect no serious degradation 
of the output of the envelope detector, so long as the time-bandwidth- 
product uncertainty relationship (equation 6-2) is not violated* 

2* OUTPUT OF THE ENVELOPE DETECTOR (SINGLE CODI; SYSTE^t) 

We v/ill investigate the single code system prior to the complementary 
coded system, since once the equations for the single code system are 
obtained, it is a simple task to obtain the equations for the complementary 
coded system* Figures 2-2a and 2-2b will be used as the system block 
diagram, since our receiver is a narrow-band system. The received signal, 
s*(t), will be given by equations 5-lOa and 5-lOb. The replica signal, 
s^|(t), is simply the transmitted signal translated by the doppler frequency; 
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hence it is identical to the narrow-band received signal. Then s.j(t) 
is given by equation 2-4. Thus 



= D'Ko'-t^ cos ( 10^4- , 



C-la 



= J(l') COS ( iOe+ t 



C-lb 



Then the output of the mismatched filter in Figure 2-2b is given by 






oo 



= tos \ Jtt-T:) (U 






J ,oo 



cos[a(i0,+uid:)1:-(u)^+t6;)t](i-b . 



-oO 



C-2 



The second term in equation C-2 will be rejected by the integrator. 
Then the mismatched output of the ‘‘envelope detector” will be given by 






f ^ 



Hli • lOA 



C-3 



From equation 6-17 we see that the mismatched output from the "envelope 
detector" is related to the wide-band ambiguity function by 



‘ 
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Recall from equation 6-14 that 

= DVD^. 

Then the differential doppler factor, in this case, is 

C-5 

Next, let us obtain the mismatched output of the "envelope detector" 
for the case of the replica signal’s doppler frequency matched to the 
received signal’s doppler frequency. This can be obtained by setting 
= 1, and 0, in equations 6-20 and 6-22. 

!• For range deviations, T , in tlie interval 






r 



Tn 






m ~ 1,2,***, 



we have 



21 

X^ir.o) = Db( i ZXiXu. 

% 

^ I 

2 [-c-v tT- 2 1 

i=X n, 1 



where N = length of the code 

X£ = itl depending on the ith value of the code symbol 
D' = received signal's doppler factor 






C-7 



= 1 + 






D. = D' 
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T == transmitted signal's duration of a code symbol 



Tj)* = received signal's duration of a code symbol 
= T/D' 

AT = magnitude of the difference in the duration of a code 
symbol = |T - 

p = smallest integer >\ ^/t\ * 



^ C-7 

(cont • ) 



J 



2, For range deviations, ^ , in the interval 



^ (m-l) IL 4 IL 4 . 

at At At at ^ 



m = 1,2, 



We have 




+2 X;Xi 

1 = 1 



Uk*i 



[r+\T- (i+u)V] 

% 







% 



M-4 
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where the definitions in equation C-7 apply and 



n = Ty'/AT 
= T/AT 

q = smallest integer 2 
k = smallest integer > 



J 



3. OUTPUT OF THE ENVELOPE DETECTOR (COMPLEMENTARY CODED SYSTEM) 

As indicated in Chapter 4-C, it is easy to obtain the equations for 
the coinpleTnentary coded system once the equations for the single code 
system are obtained by simply replacing the signed products of the single 
codes by the signed product for the complementary codes* Hence if 



in equations C-6 and C-8; then the resulting equations will be the mis- 
matched output of the envelope detector for the complementary coded 
system* 



X^Xj is replaced by XA^XAj + XE^XBj 
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AI^PENDIX D 



DURATION OF A CODE SYIIBOL DERIVED FROM CONSIDERATIONS 
OF THE SPECIAL THEORY OF RELATIVITY 

In Chapter V it was proved that the transmitted signal's duration 
of a code symbol, T, differed from the received signal's duration of a 
code symbol, T', by the inverse of the doppler factor, D. This result 
was obtained from considerations of the doppler effect on the carrier 
frequency and also by the use of the Fourier transform. 

In this Appendix we will prove the same result, but from consider- 
ations of the special theory of relativity. 

Consider the reference frame, S, v;ith the transmitter located at 
the origin of the coordinate system. Next, consider a moving frame. S', 
with the target located at the origin, A suitable set of axes is chosen 
so that the x direction of both S and S' is collinear with the velocity 
vector, Vj., At some initial time, i,e,, t = 0, x^7e shall assume that the 
two frames are coincident. At t = 0, the transmitter emits a plane wave 
pulse of duration, T. The pulse travels at the speed of light. Then we 
can view the pulse as extending in space so that its trailing edge has 
an X coordinate of -cT. For t > 0, the moving frame travels in the +x 

direction with a velocity v^. Note that the convention for positive 
velocity is away from the transmitter and is a standard convention for 
most physics textbooks. Although this convention is contrary to the one 
chosen in the main text of this thesis, no confusion should arise here. 

At time t = T, the trailing edge of the pulse reaches the origin 
of S, and the leading edge of the pulse is now at x = cT. Now let us 
consider the spacial distribution of the pulse in the moving frame. In 
the moving frame the x coordinate, x' , and the time in this frame, t'. 
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can be derived from their counterparts in the reference frame from the 



Lorentz transformations 



16 



Thus 



where 



x' = V(x - Vj.t) 
t* = 



D-l 

D-2 



Vf \ 2 T ~l/2 



i-b-ev] 



Let us denote, as in Figure D-l, the trailing edge of the pulse as event 1. 
For event 1, the x coordinate of the trailing edge of the pulse in the 
moving frame is given by 



= 



^(x - Vj.t) 



^ Vj.T . 



o 

T 



D-3 



For event 1, the time elapsed in the moving frame is given by 



l$(t --"f x) 






o 

T 



D-4 



Since the pulse is traveling at the speed of light, the total time for 
the trailing edge of the pulse to pass the observer at the origin of the 
moving frame is given by 



^ c —r~ 
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Reference frames for the 
transmitter and the target 

Figure D-1 




Reference frames for the 
target and the receiver 

Figure D-2 
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1+ Vr/c 



1 - ^r/c 



* 



D-5 



Equation D-5 is a statement of the expansion or compression of the 
duration of a pulse for the case of one-way doppler. 

Now we would like to consider the duration of the pulse which is 
reflected from the target, as seen by an observer (the receiver) at the 
origin of the reference frame, S. To avoid confusion let us denote this 
reference frame as S*'; however, bear in mind that we are assuming that 
the transmitter and receiver are at the same physical location. We can 
view the reflected wave as being retransmitted from the target. 

Consider S’, the target’s frame, and S’’, the receiver’s frame, to 
be coincident at t’ = 0. From S’, the receiver seems to be moving at a 
velocity Vj. in the negative x direction. Then at t’ =0, the leading 
edge of the pulse is at tlie origin of both frames. At t’ = T’, the 
trailing edge of the pulse is at the origin of the frame, S’. Let us 
denote this as event 2. To obtain the x coordinate of the trailing edge 
of the pulse as seen by the observer on S’’, we utilize the Lorentz trans- 
formations as given by equations D-1 and D-2, but remembering to replace 
Vj. by -v^. Thus for event 2 







k'=o 






D-6 
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For event 2 the time elapsed in S’’ 



IS 






K-O 



- ^t' 



Then the total time required for the trailing edge of the pulse to reach 
the observer on S’ ’ is given by 

T"= 




Then if equation D-5 is substituted into equation D-8 



where 



T"' 



= (^L±^k\-r 
^ 1 - ^ 

= T/D , 



0 



1- Vc 
14^r/c 



D-9 
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The factor i/hich compresses or expands the duration of the pulse is 
identical to the one obtained by the Fourier transform analysis done in 
Chapter V, with the exception of the difference in sign of the velocity 
due to the change in sign convention adapted in this Appendix, In our 
initial assumptions, we assumed that the duration of the transmitted 
pulse was T. It is easy to see that if a series of N pulses, each of 
duration T, were transmitted, then the compression or expansion factor 
would be identical to the one in equation D-9. Hence we have proved that 
the doppler effect alters the duration of a code symbol by the factor D. 
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APPENDIX E 



DETERMINATION OF THE LIMITS OF INTEGRATION 
FOR THE WIDEBAND SINGLE CODE A.MBIGUITY FUNCTION 

In this Appendix we would like to determine the limits of integration 
and their corresponding code values in the interval of integration for 
the wideband single code ambiguity function. In the formulation of the 
equations for the wideband single code ambiguity function, there are 
four quadrants to consider o We will consider only the first quadrant, 
i.e., ^>0, and will limit the validity of the equations to 

the case of the received signal’s total pulse duration being expanded 
or compressed by not more than the duration of one code symbol of the 
replica signal. 

By equation 6-17, the wideband single code ambiguity function is 
defined as 



4> i 



E-1 



Hence to determine the limits of integration and the code values in the 
interval of integration, we must analyze I(D^ t) and I(t - ) • 1(1^^ t) 

is the envelope of the received signal (time referenced to the replica 
signal), and is given by 

- O/Z ]] . E-2 



Then define 

T' = Tj)'/D . 



E-3 
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Therefore, equation E-2 becomes 






E-3 



Tlie envelope of the replica signal is given by 




E-4 



Figure E-1 is a set of waveforms of the received and replica signals’ 



envelopes 6 We have chosen, in Figure E-1, a code of length 6, with the 



we will obtain an aquation that will apply, in general, to codes of any 
length. Note that the amplitude of the received signal’s envelope in 
Figure E-1 is less than the amplitude of the replica signal’s envelope. 
This is done deliberately so that both waveforms may be seen in detail; 
however, bear in mind that both envelopes are, in actuality, of the same 
amplitude. 

For positive doppler we recall that the received signal’s waveform 
will be compressed. Figure E-1 is a series of waveforms of I(D^ t) and 
I(t), in which we hold t) fixed and shift I(t - T ) to the right, 

corresponding to positive range deviation, T . In doing so, we see that 
the limits of integration (as a general rule) change for T in increments 
of AT, where 



Next, we would like to express the interval of integration (and hence 
the limits of integration) in the form of a series. Let us now look for 



code values of 




example. From this example 



At = boi - T'\ . 



E-5 
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cont 
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•C*ToC I 1 T+4Toi 



symmetries in Figure E-1. Mote that for 04 ^ < 4 ^T, the limits of 
integration over each interval, in addition to the code values of the 
received and replica signals, remain constant. VJe also see that this 
occurs again for 9 AT4 ^ ^13 AT. However, we have lost the duration 
of one code symbol of the replica signal. Similarly, for 19AT4TA 23 AT, 
we see that we have lost the duration of two code symbols. From this 
analysis, we can obtain a general formula to express the limits of inte- 
gration and the code values in the interval of integration. Hence for T 
in the interval 






m = 1,2, •••K, 



where 



N = length of the code 

il depending on the ith value of the code symbol 
replica signal’s doppler factor 
= received signal’s doppler factor 
= differential doppler factor = D’/Di 
= T/D^ 

= 

= \TDi-T’l , 






D. = 
1 

D' = 



^Di 

T» 

AT 



r 



E-6 



The wideband single code ambiguity function becomes 


















E-7 
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in ^^7hich the definitions in equation E-6 apply and 



p = smallest integer ^ \ \ 

0( = normalization constant = ^-2^/MT 









E-8 



To complete the analysis, we would like to obtain the code values 
and limits of integration for the intervals, 4 < SAT, 14AT4^<18AT, 

and so on. As a general rule, the interval of range deviation becomes 



At at at 



wx! 

AT 



m 



“1,2,***,N* 



We see that we may express the ambiguity function for this interval as 
the sum of five terms, ioCo, 
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i-rvt 



U+Z')T' 
\ ® 





r4(i-u)To^ 



E-9 



where the definitions in equations E-6 and E-8 apply and 



n 



IL 

AT 



n 



M. 

AT 



q = smallest integer 
1 = q modulo (n ^ ) 

k = smallest integer > 






^ E-IO 
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APPENDIX F 



PPvOPERTIES OF COMPLEIIENTAPvY SEQUENCES 3, 15 

This Appendix is a compilation of some of the properties oif comple*^ 
inentary sequences » At the end of this appendix is a list of the known 
kernels, of which there are four* 

1. Definition: A pair of binary sequences of equal length, with 

the number of like pairs of one sequence, A, equal to the number of 
unlike pairs of the other sequence, B, for each possible spacing except 
at zero spacing, \;here the pairs of both sequences are like, is said to 
be a complementary sequence ^ 

2e Given sequences A and B each of length n. Then complementing 
the code, i*e«, complementing each element in the code, of A or B or 
both, results in a complementary code. 

3. Definition: Time inversing of a code is to interchange tlie 1st 

and last bit, the second and next to the last bit, and so on. 

4. Time inversing the A code, B code or both results in a comple- 
mentary sequence, 

5. Definition: The altering of a code is a transformation where 

every other bit of both codes is complemented. 

6. The result of altering a complementary sequence is again a 
complementary sequence o 

7c Definition: A kernel is a basic length code which cannot be 

decomposed into shorter length codes by an inversion of the standard 
generating methods. 

8. Definition: Complementary sequences which are not kernels are 

called composite complementary sequences. 



121 



9 • If A |Oee«eoc| 

is a complementary sequence pair, then 

C a2^,«oeoee«,a^,bj_,oe6oeoe, b^ 

B ~ a ,»e»ceeo, ,ooo*oo*, bj^ 

form a complementary sequence pair. 



lOo If A “ aj^,e«eeoe«,a^ 

B ~ b ,6«,eeeo, b^ 

is a complementary sequence pair, then 

C “ , b^ , 3 ^ 2 1 ^2 • ^n 

D — 3 -^ , bj^ , 32 9 ^ 2 1®*®*®®®®**® • ^n 

form a complementary sequence pair, 

lie The following is a list of the knovm kernels: 

a, length 2 
A = { 

B = { -1,+li 

b. length 10 

= {-l,+l,+l,-l,+l,-l,+l,+l,+l,-l} 

Bi = 'C-1,+1,+1,+1,+1,+1,+1,-1,-1,+1}- 

A2 = {+l,-l,+l,-l,+l,+l,+l,+l,-l,-l} 

B 2 = -[+l,+l,+l,+l,-l,+l,+l,-l,-l,+l} 

c# length 26 

A = -[+1,-1,+1,+1,-1,-1,+1, -1,-1, -1,-1, +1,-1, 
+1,-1,-1,-1,-1,+1,+1,-1,-1,-1,+1,-1,+1> 

B = {-l,+l,-l,-l,+i,+i,_i,+i,+i,+i,+l^_l^.l^ 
-1,-1,-1,-1,-1,+1,+1,-1,-1,-1,+1,-1,+1> 



'‘"Fi indicates the complement of 
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